Second and higher-order perturbations of a spherical spacetime 



David Brizucla. Jose M. Martm-Garci'a. and Guillcrmo A. Mcna Marugan 
Instituto de Estructura de la Materia, CSIC, Serrano 121-123, 28006 Madrid, Spain 

(Dated: February 7, 2008) 

The Gerlach and Sengupta (GS) formalism of coordinate-invariant, first-order, spherical and non- 
spherical perturbations around an arbitrary spherical spacetime is generalized to higher orders, 
focusing on second-order perturbation theory. The GS harmonics are generalized to an arbitrary 
number of indices on the unit sphere and a formula is given for their products. The formalism is 
optimized for its implementation in a computer algebra system, something that becomes essential 
in practice given the size and complexity of the equations. All evolution equations for the second- 
order perturbations, as well as the conservation equations for the energy-momentum tensor at this 
perturbation order, are given in covariant form, in Regge- Wheeler gauge. 
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I. INTRODUCTION 

When a physical problem cannot be solved exactly, one 
usually appeals to approximate methods. Perturbation 
theory provides one such approximate approach, allow- 
ing a description in terms of small departures around an 
exact solution. In the context of General Relativity, per- 
turbation theory plays a prominent role in analyzing and 
understanding dynamical processes, being nowadays an 
efficient and natural complement to full numerical rela- 
tivity simulations 0. 

In particular, perturbation theory is used to study the 
stability properties of solutions of interest: black hole 
spacetimes cosmological solutions 0, critical solu- 
tions 0, and many other. It also allows us to check the 
resence of gauge instabilities 0, constraint violations 
, and other types of instabilities in the various formula- 
tions of the Einstein equations implemented in numerical 
relativity, since numerical errors can be considered them- 
selves as distortions of the solution that one is computing. 
Of outmost importance, perturbation theory can provide 
us with estimates of the amount of gravitational radia- 
tion and of the signal profiles emitted in astrophysical 
scenarios like an oscillating neutron star 0, the gravita- 
tional collapse of a star 0| , an extreme mass-ratio binary 
13, or a close limit head-on collision of two black holes 

The complexity of the expressions involved in the per- 
turbation of the equations of General Relativity is very 
high, rapidly increasing when working at higher orders of 
perturbation. Most of the previous investigations have 
been carried out at first order, and using highly sym- 
metric unperturbed solutions, what simplifies the prob- 
lem. Simple backgrounds have been perturbed at second- 
order level: for example, the Friedmann-Robertson- 
Walker spacetime the Kerr spacetime O'" the 
Schwarzschild spacetime 0,^3 • The change in the oscil- 
lation modes of a stationary star owing to its rotation has 
been studied in Refs. regarding the rotation as a 

first-order perturbation of a static star (and therefore the 
change as a second-order effect). The critical exponent of 
angular momentum scaling has been predicted for scalar 



field collapse using second-order perturbation arguments 
|l7j . Special credit deserves the seminal work of Cun- 
ningham, Price, and Moncrief in 1980 where the second- 
order nonspherical perturbations of a collapsing star of 
dust were studied including the problem of matching of 
the internal and external perturbations through the sur- 
face of the star ^3 • Second-order perturbations of the 
general problem of matching through a surface have been 
analyzed in Ref. [T^ . 

Different reasons justify the importance of going be- 
yond first-order perturbation theory. First, one must 
confront the obvious desire to reach higher accuracy in 
the numerical simulations of perturbative approxima- 
tions to self-gravitating systems. In addition, second- 
order perturbations can estimate quantitatively the range 
of applicability of the first-order results by providing er- 
ror estimates. Besides, they should enable us to study 
and interpret the nonlinearity of General Relati vity in 
terms of the coupling among first-order modes [23|, a 
study that might allow us to model secular interactions 
which arc too slow to be followed by using full numerical 
simulations or which might be incorrectly interpreted as 
small numerical errors. 

In this work, we construct a generic framework to ana- 
lyze second and higher-order spherical and nonspherical 
perturbations of an arbitrary spherical spacetime, with- 
out restricting to any particular matter model, and show 
how to deal with the general case while keeping a still 
manageable size for the resulting expressions. This work 
can be considered as a continuation of the work of Ger- 
lach and Sengupta (GS) at first order a for- 
malism which was revived in Ref. |23l | and is currently 
considered optimal for perturbations of generic spherical 
spacetimes [24Ll25j (even though for specific matter mod- 
els the formalism might admit further simplifications). 
For example, its application to a perfect fluid spacetime 
with a two-parameters equation of state was considered 
in Ref. , including the matching of the fluid perturba- 
tions to an exterior spacetime through a moving timelike 
surface |2^ . 

The GS formalism is based on four basic ingredi- 
ents: a 2+2 decomposition of the spacetime separat- 
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ing the spherical S'^ symmetry orbits from a general 
1+1 Lorentzian manifold Af^; the use of a covariant de- 
scription both on the Lorentzian manifold and on the 
2-sphere; the decomposition of the perturbations in 5*^ 
tensor harmonics; and the use of gauge-invariant pertur- 
bation variables. The formalism developed in the present 
paper makes use of the first three ingredients, leaving 
the construction of gauge-invariants for a future work. 
The use of a covariant notation is particularly convenient: 
on the one hand, it allows us to formulate all equations 
without choosing coordinates on M^, something that be- 
comes very useful on dynamical backgrounds; on the 
other hand, it eliminates all trigonometric factors from 
the equations of motion, factors which do not contain any 
relevant information and typically obscure the geometri- 
cal interpretation of the results. The use of covariant no- 
tation, however, makes computations more complicated. 
To overcome this problem we will intensively use com- 
puter algebra tools, specially designed and developed by 
us for abstract tensor computations. All equations in this 
paper have been programmed and checked with those 
tools, resulting in a computer framework which allows us 
to work efficiently with applications of the formalism of 
high-order perturbations. 

For some matter contents it is possible to further sim- 
plify the perturbative formalism by constructing scalar 
combinations of the perturbations which encode the 
purely dynamical degrees of freedom and obey evolu- 
tion equations free of constraints (the so-called master 
equations for the master scalars). Among other matter 
models, this is possible for vacuum [28L |2^ and for the 
Maxwell field iSQl, but such master equations have not 
been found, for instance, for the case of a scalar field. 
Those scalars can be given in a gauge-invariant form [sij 
and expressed in covariant form within the GS framework 
[21L [33 . In this article we will concentrate on the con- 
struction of a generic perturbative formalism that can be 
applied to any background spherical spacetime. The con- 
struction of gauge-invariant master scalars at high per- 
turbative orders will be the subject of a future work. 

The rest of this article is organized as follows. Section 
II introduces the fundamental concepts and equations of 
perturbation theory, and provides closed formulas for the 
nth-order perturbations of the geometric quantities of in- 
terest (formulas that are new in the literature to the best 
of our knowledge). The notation used for the spheri- 
cal background manifold is explained in Section III. Sec- 
tion IV introduces the Regge- Wheeler- Zerilli harmonics 
|28l Ib^ I and generalize them to arbitrary rank tensors. 
Formulas for their products are given, based on the repre- 
sentation matrices of the rotation group. Basically, noth- 
ing in that section is new, but it has been conveniently 
recast into GS notation, which proves to be very useful 
for applications and computational efficiency. Nonspher- 
ical perturbations are discussed in Section V, giving for 
the first time complete sources for the evolution equa- 
tions of second-order perturbations of a general spherical 
background, as well as sources for the energy-momentum 



conservation equations. Section VI contains a summary 
of our results and further discussions. Finally, several 
appendices are added. They explain different aspects 
of the definition of spherical functions, the symmetric 
trace- free part of tensors, and the pure-orbital harmon- 
ics. They also provide the GS equations for first-order 
perturbation theory and describe the procedure followed 
to implement our calculations in Mathematica [3^. 

II. PERTURBATION THEORY IN GENERAL 
RELATIVITY 

In General Relativity, one has to solve ten coupled non- 
linear partial differential equations for the metric. Using 
perturbation theory, the problem can be reformulated as 
an infinite hierarchy of linear differential equations for 
the modification of the metric with respect to a known 
solution. Furthermore, in this hierarchy of equations, the 
principal parts are always given by the same differential 
operator acting on a perturbative correction of increasing 
order. In this way, one translates the difficulty from non- 
linearity to the infinite number of equations. The key 
assumption of the perturbative scheme is that one can 
truncate the problem at a finite order and still obtain an 
approximate solution to the original system. 

In order to introduce this perturbative hierarchy, let us 
start by considering an e-family of 4-diniensional metrics 
Qiiui^) on a certain manifold Al, where e is a dimension- 
less parameter jZ^I- The matter field content of these 
spacetimes provides another e-family which will be de- 
noted abstractly by ^(e). For simplicity, we assume that 
these families are smooth in the parameter e (or at least 
C" with a locally Lipschitz n-th derivative if we are inter- 
ested only in perturbation theory up to order n > 0). In 
particular, the metric and matter fields can be expanded 
as 

Qt^A^) = 9^.1' + ^ — ^''^''h^u, (1) 

n=l 
00 „ 

$(6) = ci> + ^i^w$. (2) 
n! 

n=l 

The e = fields and $ will be referred to as the 
"background" metric and matter. To distinguish the 
background objects from their "perturbed" e 7^ coun- 
terparts, we will denote the latter with a tilde (""). 

The coefficients and are tensors on the 

manifold Ai and in what follows their indices are 
lowered or raised with the background metric g^j^i, and 
its inverse g^''. They can be obtained by repeated dif- 
ferentiation of the perturbed objects with respect to e. 
Actually, it proves most convenient to introduce a for- 
mal "perturbation" operator A with the properties of a 
derivative, so that any object T(e) can be expanded as 

00 „ 

fie)=T + J2'-A-[T]. (3) 

n— 1 
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For instance, A [g^,^] = and A[ = + 

In this notation, the brackets are intended to avoid con- 
fusion with index positioning because, e.g., for a given 
vector u^: 

(4) 

so that the notation Aw^ might be misleading. 

The expansion for the inverse metric can be obtained 
by iteration of the identity 

^g'-'' -a^-^isxa-axaYa"" (5) 

or, equivalently, by repeated perturbation of the relation 
A[5^1 - -5^"A[g„^]/^ (6) 

This leads to 

and, more generally, to a perturbation of the form 

(8) 



fci! ... fc,„! 



where the sum extends to the 2"~^ sorted partitions of 
n in m < n positive integers fci + ... + km ~ n. For 
example, for n — A there are eight partitions: (4), (1,3), 
(3,1), (2,2), (1,1,2), (1,2,1), (2,1,1), and (1,1,1,1). 
The Christoffel symbols are 



(9) 



where we have defined the three-indices perturbation 
{"}/, =i (-{"}/, -I- (">/; -<">/) ^ fin) 

which is symmetric in its last two indices and satisfies 
^"^hafii, = 0. The covariant derivative in H1Q|) is that as- 
sociated with the background metric. Higher-order terms 
of the expansion can be easily computed noting that 

A [ = f" + ^>/i„^, - '"'hj (11) 

which, for n > 1, leads to 

nl 



A"[r"^,] = ^(-1)"+^ 



. ... h,Yn • 



(12) 



Here, the sum extends again to all sorted partitions of n. 
The case n = 1 is special, with A[r"^y] = Note 
that each term in the above expression contains one and 



only one tensor ^''^/ip^i/, but a variable number of metric 
perturbations. 

The perturbations of the Riemann tensor are given by 



Rn 



— R^iua ^" 2e ^ (t[/^;i/] 



(13) 



+2<^>/i„[/f^>/i%],)+0(e3), 
with the general term 

A"[i?p,„^] = (14) 



fc=l 



(m ^ v). 



This expression does not involve the metric directly. 
That is, it only contains the background connection V, 
and the perturbations of its Christoffel symbols, without 
assuming that either the background or the perturbed 
connections derive from a metric. Hence, it can be ap- 
plied to the Palatini equations, for example. If the back- 
ground connection derives from a metric, Eqs. (jSl^^l 
ensure that the perturbed connection also derives from a 
metric. Then, the (implicit) triple sum in Eq. p4fl can 
be rearranged as follows: 



A"[i?M./] = E(-l) 



/si!...fc„J 



(15) 



rn 
s=2 

-(m ^ v)- 

This formula is surprisingly simple because all covari- 
ant derivatives of the metric perturbations are grouped 
in /i-terms of the form (I10|l . It is clear that the explicit 
sum in Eq. (|14|l contains only that kind of terms; how- 
ever the derivatives of the perturbations of the Christoffel 
symbols give rise to isolated covariant derivatives of the 
metric perturbations. Nonetheless, they can be combined 
with the sum to obtain the displayed result. Note also 
that this expression is already optimally simplified: all 
terms in the sums are generically different. 

From the above formulas, we can compute the gen- 
eral perturbation of the Ricci and Einstein tensors. For 
Ricci, the perturbation is obtained from Eq. (|15|1 by con- 
tracting the (3 and v indices. The resulting expression is 
symmetric in a and ^ owing to the identity 

n n 
^ ^h^^^,^^ = ^ WT'^P '"'Kp,,^, (16) 



k=l 



k=l 



valid for any family of symmetric tensors f'^'^T"^ such 
that X]fc=i ^''^T°'^ is also symmetric in a and ^. 
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It is important to emphasize that these combinatorial 
formulas for the perturbations of the curvature tensors 
at a general order are extremely useful for computational 
purposes. Essentially, the problem of perturbations is re- 
duced to that of listing the sorted partitions of a given 
number, which can be done very fast in any computer- 
algebra system. This fact simply reflects the recursive 
differential origin of the perturbation process. The for- 
mulas in this section have been implemented in the free 
package xPert, briefly described in Appendix lEl 



III. SPHERICAL BACKGROUND 



where eg" and e^" are the unit norm (with respect to 
the round metric) basis vectors. From these vectors we 
obtain 



a relation which can be inverted to get 

e'^f' = i{m''m^ -rh'^m''). 



(21) 



(22) 
(23) 



These vectors are null, "fabfn'^iTL^ — '~fab'm'°"m-^ = 0, and 
are normalized so that ^ab'm°"fn^ = 1- 



Following GS [22, |23|, the manifold of a spherical 
spacetime (M^, g^^) is described as a product of the form 
M'^ = X using a coordinate system a;'' = (x^, a:") 
adapted to the 5^ orbits of spherical symmetry. The 
two scalars x^ = {x^,x^} provide coordinates for the 
1+1 Lorentzian manifold with boundary, whereas 



{x^ 



i,X-' 



4>\ are the usual spherical coordinates 



on the sphere S'^ . The 4-dimensional metric g^^i, and the 
energy-momentum tensor t^^, can always be written 

g^,^{x^,x'^)dx''dx'' = (17) 

gAB{x'^)dx'^dx^ + r^{x^)-fab(.x'^)dx''dx\ 
tf,^{x",x'^)dx^'dx'' = (18) 

tAB{x^)dx^dx^ + K^{x^)Qix^)-fab{x'^)dx''dx\ 

where r is a scalar field on , gAB is a Lorentzian met- 
ric tensor on M^, and jab is the round (unit Gaussian 
curvature) metric on S^. Using this decomposition it 
is straightforward to express all 4-dimensional curvature 
tensors in terms of the curvature tensors of g and 7 and 
the derivatives of the scalar r, as shown in Ref. |23|. In 
doing so, one usually introduces the vector field 



VA = = (fog7'),A 

r 



(19) 



to avoid working with logarithms of r. Covariant deriva- 
tives of a vector / on the manifolds M^, AP and will 
be denoted as fA\B and fa-.b, respectively. For in- 

stance, the Schwarzschild spacetime has been described 
along these lines in Refs. [Mill HI. 

The totally antisymmetric tensors are denoted with 
the symbol e and obey the conventions CABcd = ^AB^cd 
with £01 = 1 on and £23 = 1 on 5^. Note that 
some authors interchange the uppercase/lowercase index 
conventions, and others use the opposite sign for eqi. 

On the other hand, as remarked by Newman and Pen- 
rose [35I I , it is convenient to introduce a basis of complex 
vectors on S*^: 



V2 



(ee" -f ie^"), and m° 



V2 



IV. TENSOR HARMONICS 

The theory of tensor harmonics on the sphere S*^ is 
a very well known subject js^ Is^, 113 ■ In this section 
we recast it into the GS notation, which is particularly 
useful for the type of algebraic computations that must 
be performed in high-order perturbation theory. After a 
brief review of the GS harmonics we show that, in or- 
der to handle second and higher orders of perturbation, 
it is convenient to work with harmonics that posses an 
increasingly high number of indices. A generalization of 
the GS harmonics to many indices is then defined and 
shown to be closely related to the Wigner rotation ma- 
trices (also known as spin-weighted harmonics in General 
Relativity js^ ) . Finally, we give closed formulas for arbi- 
trary products of these harmonics. Appendix fXl explains 
our choice of conventions. 



GS notation for harmonics 



An orthonormal basis of functions on the sphere S"^ 
is given by the spherical harmonics yj™(6', </>), which are 
defined as the eigenfunctions of the differential operators 



r'^^a^bY^-- 



-i{i + i)r/' 



(24) 
(25) 



where 7"^ is the inverse of the round metric on S'^, the 
angle cf) is defined by choosing a fixed z-axis, and I and 
TO are integers such that I > \m\. These harmonics are 
normalized so that 



dnr,:' 



Si'i s„ 



(26) 



(20) 



with dil being the area element on {d^l = sin 9 dO d4)) . 

From them, the Regge- Wheeler (RW) basis of vec- 
tor fields on the sphere can be defined as follows. The 
basis is formed by the vector fields Y^^-a and their or- 
thogonal fields Sl^a = ^abJ^'^Y/^.c- A basis for 2-tensor 
fields on the sphere can also be constructed in a simi- 
lar way and it is formed by three types of objects: 
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pure-trace tensors can be decomposed using "fabYi"^', an- 
tisymmetric tensors can be expanded using Cafc^™; and 
finally symmetric traceless tensors can be expanded using 



; ab 
ab 



(yr-.ab) 

(a:b) 



TF 



Y, 



I -ab 



labYr, (27) 
(28) 



Indices in round brackets are symmetrized, and the su- 
perscript TF denotes the trace- free part. Note that GS 
use 2Sl^i^a:b) instead of XPab, different by a factor of 2. 



B. Products of harmonics 

In the next section wc will expand the metric perturba- 
tions '"'/i^^ in tensor harmonics. From expressions like 
(jHl; H12|l . and 1)15(1 it is clear that wc need to compute 
products of several tensor harmonics when working be- 
yond linear perturbation theory. Even though those ex- 
pressions contain products of many harmonics, the prob- 
lem can be dealt with recursively because the product 
of two tensor harmonics can be decomposed as a series 
of tensor harmonics of adequate rank. In principle, we 
might conclude that at perturbation order n we need to 
work with tensor harmonics of rank 2n or similar, but the 
situation turns out to be simpler in General Relativity. 

The formalism starts from perturbations of the metric, 
which contain tensor harmonics on S*^ of rank zero, one 
or two, and computes the decomposition in harmonics of 
the perturbations of the Einstein tensor, which also con- 
tain harmonics of those ranks. On the other hand, only 
second (at most) derivatives of the metric perturbations 
will appear in the perturbations of any curvature tensor, 
at any order. Finally, as long as we are interested just 
in perturbations of curvature tensors, only those contrac- 
tions in Eq. H15() are required. From these three observa- 
tions wc conclude that we only need harmonics with up 
to four indices (if one works in RW gauge, to be defined 
below, only three-index harmonics are required) and for- 
mulas for their thirteen products 













YY'..a, 


YaY' 


bi 






YY'.,ab, 


YaY' 


be •} 


Y..abY' 


■.cdi 


YY'..abc, 


YaY' 


bed 


YabY' 


■.bed-) 


-.abed-) 


YaY' 


bede : 


YabY' 


■.cdef 



Y-.abcY -.deft 

(29) 

where the prime denotes that Y and Y' have different la- 
bels I and m. Only seven of those are really independent 
because using the Leibnitz rule we have relations like 



Y..abY'..ed = {YaY'..cd):b - YaY'., 



cdb ■ 



(30) 



Therefore, computing the expansion formula for the 
canonical products YY' ■aj^...a„ with n = 0, ...,6 would be 
enough to solve a general problem of nonspherical per- 
turbations in General Relativity. 

That method would be, however, rather complicated to 
program, because it requires to expand the products of 



multiple harmonics in a very particular order, and diffi- 
cult to use in any mathematical proof involving products 
of harmonics. It is far more interesting and general to 
follow a different route: we first generalize the GS har- 
monics to an arbitrary number of indices and then find a 
general formula for the product of any two of them. This 
has two important advantages: first, it is more efficient 
and simple for our algebraic code because all cases are 
considered in a single formula. Second, the formalism is 
more general: it can be applied to arbitrary matter mod- 
els, it is possible to perturb objects like derivatives of the 
Riemann tensor, or can be used in other problems (for ex- 
ample theories of gravity with more than two derivatives 
in their basic equations). 



C. Higher-order generaHzation of GS tensors 

Complete bases for second and higher-order tensors 
can be easily constructed. There always exist two non- 
trivial symmetric trace-free (STF) tensors 



{Yi"^:ai...aJ 



STF 



ai...as = W/ ai:a2...asj 



\STF 



-^{ai^XPba2...as},{'^^) 
■ e(ai^^rha2...a,),(32) 



valid for \m\ < I and 1 < s < I. In all other cases the 
harmonics arc defined to be identically zero, except for 
s — 0, when = YJ™. Note that, in fact, we do not 
need symmetrization on the far right-hand side because 
the tensors Z and X are traceless. All other objects in 
the basis can be obtained from products of 7, e and the 
basis for tensors of order s — 2. For example the basis 
for 3-index tensors is given by ZJ^abc, X^abc, and six 
independent combinations of ^abZf^e, JabX^c, ^abZ^c, 
eabXl^c, and their index-permutations. The general case 
results from the iteration of the relations (valid for s > 2): 



/ ai...as:b — ai...asb v 



{l + s){l-s + l) 



i{a\a2^l a3...as)b fb{ai'^l a2...as) 



(33) 



vm I 

ai...as:b — ai...asb i 

'1 



(34) 



-,l{aia2^l a3...a^)b — lb{ax^l a2...a^) 



Appendix Ogives a different approach to expand the def- 
initions lEU and (|^ . 

Remembering the definitions of the scalars = Y[^ 
and X™ = 0, and those of the vectors Z]^a = Y^^.a and 
Ar™a = 5'["a, we obtain the three remaining special cases: 



I -.a — a, 



7m yin 

-H a:b — ab 



x^ 



I a:b — ab 



1(1 + 1) 
2 

1(1 + 1) 





(35) 


'7™ 
labZi , 


(36) 


-^ab^l ■ 


(37) 
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Formulas |(221l~lE7I) for the STF tensors Z and X con- 
stitute a complete set of simplification rules which allow 
us to express any derivative of a tensor harmonic field 
on the sphere in a unique canonical way. Note that ^ah 
appears only in Eq. H37() . 

Finally, it is important to point out that all harmonics 
have a well-defined parity under inversion of axes. This 
is because the parity of 7af,, tah^ and the scalar harmon- 
ics Y;™ is -t-1, -1 and (—1)', respectively, and because 
taking covariant derivatives does not change the parity. 
It is usual to separate harmonics with momentum / into 



two families: tensors ZV 



have parity (—1) (these 



are said to be "polar" or of "even" polarity) , and tensors 
XPai...a, have parity (—1)'+^ ("axial" or "odd" polar- 
ity). One must not confuse parity and polarity: whereas 
all equations must have a well-defined parity at any or- 
der in perturbation theory, polarity is only useful in the 
first-order theory because products of harmonics couple 
the two polarities, as we will see. 



D. A formula for the product of tensor harmonics 

The product of two scalar harmonics can be expanded 
in terms of finite sums of scalar harmonics using Clebsch- 
Gordan coefficients |3^ : 

I'+i 

T^m' v^m \ ^ 77^0/ m \r7n+m' fOQ\ 

l" = \l'-t\ 

where we have defined the symbol 



E 



'01 m 
Ol'm'l' 



_ / (2? + f )(2?^ + 1) ^m'mm'+m/^OOO 

= V 4^(2?" + 1) ' 



(39) 



Recall that the Clebsch-Gordan coefficient C°"°/ van- 
ishes if r -I- ^ -I- 1" is odd. This fact guarantees that only 
scalars with parity (—1)' = (—1)' are present in the 
expansion. 

In this subsection we will construct a generalization of 
Eq. (|38|l valid for any pair of tensor harmonics on the 2- 
sphere. There are two main routes to find such a formula. 
The standard route, explored in Appendix|21and followed 
by most books in Quantum Mechanics, is adapted to the 
3-dimcnsional Euclidean structure of and uses the so- 
called "pure-orbital" harmonics C'/'™ii...i3 , which trans- 
form under a representation of "total angular momen- 
tum" J, with \m\ < j, and whose Cartesian components 
are eigenfunctions of the "orbital angular momentum" 
operator l|C3|) with eigenvalue l{l + 1). This latter prop- 
erty becomes very useful when solving wave equations in 
a 3-dimcnsional setting. Unfortunately, these harmonics 
are not transverse to the radial direction, a fact that un- 
necessarily complicates the analysis of the radiation in 
the far region. 

Here we will follow the second route, based on the 



so-called "pure-spin" harmonics 3^, 



±s,m 



They are 



adapted to the 2-sphere, and hence are transverse to the 
radial direction. Besides, they are closely related to the 
Wigner representation matrices of the rotation group, for 
which a product formula is well known. These harmonics 
can be defined in the following way. 

Let us consider the unit sphere, with points described 
in a certain frame by coordinates {6',(j)'). A rotation R 
of that frame assigns new coordinates {9, (f>) to the same 
physical points, so that scalar fields / transform as 



f{0,(t>)^f{e',ci,') = D{R)f{e,, 



(40) 



Since spherical harmonics provide an irreducible repre- 
sentation of the rotation group, their transformation rule 
must have the form 



DiR)Y/^ 



I 

7n' — ~l 



(41) 



where a, /3, and 7 are the three Euler angles correspond- 
ing to the rotation R. Explicit formulas for the compo- 
nents of the matrices ^'^^/^(q!, /3, 7) have been computed 
by Wigner and are reproduced in Appendix^ Most 
important for us, the product of any two of these ma- 
trices for the same rotation R can be expanded using 
Clebsch-Gordan coefficients: 



V 



]' (R)V'-^V (R) 



(42) 



TOi m2 mi +m2 ™2 ™'l +™'2Tt{i) 



32 3 



31 32 3 



V''> , ^ (R). 

771^^+7712, mi +7712 



Making use of these matrices, the spin- weighted scalar 
harmonics can be defined as 1391 



s ^ L,m : ' 



21 + 1 



(0,0,0). (43) 



In particular oYi^,n = Yi"^- 

STF tensors on the unit sphere can be constructed 
from the vectors m" and to" defined in Eq. H20() . Two 
independent STF tensors of rank s arc 

to"i...to"= and m.''\..m°-% (44) 

because ^abTn°"m^ = = ^ab'fn""ff^^ ■ We then define the 
pure-spin tensor harmonics with s > indices on as 

yt"\^...a^ = (-l)"A:(Z,s)2?(^)„(O,0,0)TOai...TOa,,(45) 

3^r''"ai...a. = k[l,s)V%^{Q,eA)rha,-rha^, (46) 
with 



fc(/,s) 



'{2l + l)[l + s)\ 
2''+2 7r(/- s)!' 



(47) 



The normalization factors arc introduced so that the GS 
harmonics arc (for s > 1) 



I ai...as 



I ai...as 



y, 



— s.7n 

I ai...as 



(48) 
(49) 
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,0,m 



For the case s = 0, one has = y[ 

We can invert the previous relations to get (except for 
the special case s = 0) 



I ai...a.. 



(mQim''y™:6Q2...aJ 



STF 



■.ba2...a^ 



, STF 



(50) 
(51) 



Formula 142(1 provides the following product of pure-spin 
harmonics with the same sign: 



must be expanded using Eq. H21|l . We define T° = 1. 

Adopting the notation <='='Z = y^zty^'^ for all s (where 
we have obviated the rest of sub and superindices), and 
introducing the tensors = ^{T'" ± T'^) and the al- 
ternating sign e = (—1)'+' ~' , we hence arrive at the 
final formula (assuming again that s' > s) valid for the 
product of any two generalized harmonics 



J^lf ai...a^rJ^i bi...bs — 

I'+l 

ETP±s'l'm' ,,±(s'+s),m'+m 
^±s I m l"-yi" ai...a^,bi...b^ 

l" = \l-l'\ 

where we have introduced the real coefficients 



(52) 



^s'l'm'l" 



k{l , |g Df^i^A^l) rtrn'mm'+mris' ss'+s 



(53) 



which generalize the coefficients (|39|l . These inherit from 
the Clebsch-Gordan coefficients the symmetry properties 

i?z:'r;;:\. - Et' fi™',' = {-iy'+'-'" Ei'\'z(,,, (54) 

(55) 



rps' I' m' 
I m /" 



^s'l'm'l" 



From the fact that CpY'f^ 



for odd I" we also get 
that the i?-coefficients vanish for odd /" ii I = I' and 
either m = m' or s = s'. 

For the remaining products of pure-spin harmonics 
(those with opposite signs), we obtain (assuming e.g. 
that s' > s without loss of generality) 



;'+/ 

/ m ai...a^rbi. 

l" = \l'-l\ 
I'+l 



(59) 



t \ ^ ^T?s' I'm' ( {Eaa')'j'm' +irt ^7"+ 
+ ^^-slmVy ^l" a.+ l---a,, ■'aibi...a,b, 

-I- ^in as + i...a^,-'aibi...asbsj ' 

which constitutes the main result of this section. The 
first sum in this formula is very simple [similar to that in 
Eq. (|52() ] and involves only harmonics with s' + s indices. 
The second sum involves harmonics with s' — s indices 
and has a more complicated structure in order to include 
the case of products with scalar harmonics. 



(56) 



~ii=ps ,?n -lilts, m 

J^j' ai...a^,.yi tii...bs " 

/'+( 

Erp^s'l'm' ->,t(s'-.'5),™'+™ 7^±s 
; m l"-yi" as + i...a^, -i- aibi...asbs 



NONSPHERICAL PERTURBATIONS 



A. 2+2 decomposition 



where the products 

a\b\...asbs 

T—s 
a\b\...asbs 



{~iymaimb^...ma^mb^, (57) 
{-lYmairhb^...'ma^mb^ (58) 



From now on, we adopt the following 2-1-2 decompo- 
sition of the perturbations of the metric g^^ and the 
energy- momentum tensor i^i^: 



I ab 



(60) 



/it/ 



E 

Lrn 



I A ^. 



I AB 



vm 
I A^l b 



r 



""labzr - 



I a.b 



I b 



I a.b 



(61) 



The polar (axial) components of the perturbations are have been displayed in this expression, forcing us to use 
denoted with uppercase (lowercase) letters. All indices 
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up to five indices in the tensor perturbation ^'^''HY^ab- 
For n = 1 these expansions in harmonics (and in partic- 
ular the choice of factors ) reduce to those of Ref. [23| , 
with some changes in the notation and the mentioned 
difference of normahzation of the axial tensor Xf^ab- 



C. Evolution equations 

The spherical background metric H6U|) satisfies the Ein- 
stein equations G^u = 87rt^i^, which can be decomposed 
as El 



B. Gauge dependence 

The perturbations (|60|l and H61|l are not invariant un- 
der changes of gauge, i.e., under changes of the point- 
to-point identification between the perturbed and unper- 
turbed spacetimes. Those changes can be thought of as 
diffeor norp hisms on one of those manifolds. Soncgo and 
Bruni |40| have given a complete description of the Tay- 
lor expansion of families of diffeomorphisms around the 
identity map. A given family is described by an infinite 
collection of vector fields <">^'' -one at each order in the 
power series expansion- in terms of which the action of 
the diffeomorphisms can be expressed using Lie deriva- 
tives. For example, up to second order the result is 



AfTl 



A[T] 
A2[f] - A^iT] 



CiD^T, (62) 
C in^T (63) 

2/:{i,.A[r] + /:2„^T, 



for any tensor field T, where the ovcrline denotes a dif- 
ferent choice of gauge. In particular, for the metric field. 



{1}/, _ {1}/, 



= C{i)^gt,y, (64) 
= C{2}^g^^ (65) 
-t- 2C {1}^ + C\n^g^^. 



These formulas allow us to compute the perturbations in 
any desired gauge from their values in a particular one. 

In principle, it may be possible to form gauge-invariant 
combinations of the perturbations at any desired order. 
For simplicity, nonetheless, here we will perform our cal- 
culations by imposing a particular choice of gauge: the 
generalization to higher orders (n > 1) of the RW gauge, 
namely. 



I A 



0, 



0. 



(66) 



This leads to a full metric 5^1/ whose components Ab and 
ab obey four local gauge conditions at all points, 



gAb:c9 = 0, 



gab = Kgab, 



(67) 



for some generic scalar field K on M^. To be more pre- 
cise, the RW conditions select the point-to-point map- 
ping between the perturbed and unperturbed spacetimes 
in such a way that the pull-back g of the perturbed met- 
ric into the background manifold obeys Eq. (|67|l . Con- 
versely, using these four conditions it is easy to see that 
it is always possible to impose the RW gauge at all per- 
turbation orders in any family of metrics, at least at a 
local level, so that the gauge is well posed. 



Gai 



Ga 
a 



'2{vA\B + vavb) 



(68) 



9AB {^-^ + 2vc^'- +3vcv'' 
-'^''R+2vA^^ +2vAv'^ = 87r( 



87r tAB, 



(69) 



The evolution of the first-order perturbations is sche- 
matically given by the six GS equations 



(70) 
(71) 
(72) 
(73) 



Eab[ 




- 8^*^>*rAB, 


Ea[ 




= Stt^'^^^a, 


E[ 






E[ 






Oa[ 




= 8^(^f^vrA 


0[ 




= 8^^^vr, 



^Q'^'hrA 



(74) 
(75) 



where represents the 10 first-order perturbations 

with labels / and to, and the E and O differential opera- 
tors are given in Appendix IdI The E operators contain 
only polar metric perturbations and the O operators only 
axial ones. 

The evolution of the second-order perturbations is dic- 
tated by the same equations, except for that now the 
left-hand side contains extra sources that are quadratic 
in the first-order perturbations: 

EABi'^'hr] <^'5pf rAB-87rf^>vI/-^B,(76) 

Li 



EAi^'^h]^] + ^ ^ '^'^pfr^ = 87rf^>*rA, (77) 

l i m,m 

Ei^'^h^] + ^'^Sf'fr = Sn^'^^^r, (78) 



I I m,7n 



Li 



mmrn 



37rf^>*[", (79) 



Oa^ 



( - e ) c m mm 



(80) 



Btt 



{2},/,™ 



■ / m 
Vl A 



^EE 



<--)S'pt"J" = 87r*=Vr, (81) 



with the usual restrictions on the values of m and m, and 
with both I and I being independent and running over 
all nonnegative integers. The structure of the sources is 
rather peculiar, owing to the mixture of polarities that 
appears in the product of harmonics. This fact is encoded 
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in the polarity sign a of the sources '"'Sppj", which is al- 
ways given in terms of the associated sign e = (—1)'+'^' 
and thus completely determined for each term of the sum, 
so that it cannot be chosen freely. Sources with polarity 
sign CT = +1 contain terms polar x polar and axial x axial 
with real coefficients. Sources with polarity sign a ~ —1 
contain terms of the form polar x axial with purely imag- 
inary coefficients. This form ensures an adequate behav- 
ior of the equations under complex conjugation (that we 
denote with the symbol *). In particular, using 1)54(1 . we 



have for all sources and for all LLl: 



I i I 



{ — £) grti rn rn 



I i I 



(e) Q~m—m — m 



) Q — vn—m—m 



(82) 
(83) 



The sign e alternates when any of the I labels changes. 
Therefore, all equations have generically both types of 
sources. 

On the other hand, we see that some pairs of equa- 
tions share the sources: for example, Eqs. 1(79(1 and l(81() 
alternate their sources '+'5* and '^'S* for particular sets 
of labels 1,1,1. The same thing happens with the pair 
(177(1 and ((80(1 . The operators Eab and E, however, have 



their own pair of sources. As a result, we need to com- 
pute eight sources in total, instead of twelve. 

Using the expansion ((15() and the definition of the met- 
ric perturbations ((60(1 . we can expand the Einstein equa- 
tions at second order assuming, without loss of gener- 
ality, that there are only two first-order perturbations, 
but allowing these to be completely arbitrary, in partic- 
ular assigning arbitrary harmonic labels to them. From 
now on the coefficients and harmonic labels of those two 
perturbations will be denoted as h and h, with all other 
perturbation amplitudes vanishing. In this way, we avoid 
to deal with sums that include (quadratic) couplings 
between an infinite number of first-order perturbations. 
The expansion contains many terms with products of ten- 
sor harmonics: we count 1275, 972 and 1347 source terms 
in A[Gab], '^[GAb] and A [Gab], respectively (still at the 
24-2 abstract level, without any expansion in coordinate 
ranges) . These products of harmonics must then be ex- 
panded using formula ((59(1 . We now analyze the sources 
separately. 

The source of A [Gab] contains products of harmonics 
of the form Z'^'"^Zabc, Z°'''Xab, etc. (harmonics with four 
indices do not appear in RW gauge) . The final expression 
can be rearranged to arrive at the sources: 



(+) nramm 

I I AB 



— ^2 irhl9ABtlcn ^ 2 5 ^0 [rhl ['T-AhB - 9AB hc fl 



(84) 



1 



2^1 ^%\^f^'^ihA\B)C - hc\AB + hcVA\B + '2hc\{AVB))+^h'^\ch(A\B) " 2h^ \Ahc\B 



2hA^^hB\c + 9ab 
2 '^'i? - 4 



AvcVD]-2r-\r'h^)\cir'h^)\n 



r^ 



r^ 



1 



2,-2 1 f™' 



2H^'cHab - ^HacH^b + 9AB UkK + SH^'^Hcd - H^cH^d 



P+l + P+l-2 



HabK + H^abK\c + -HABir^Kic)^^ + {KK)^ab - K\aK\b 



1 

2' 

tF 



C 



Av(aK\b)K + -Hcd(aH^b}'' - H^"" [r-Hr^HcAB)iD - Hcd\ba] - -HdabH^^'c 

H^c^^D - H^^icD + 2H^c\DV^ - AH^^icVD - 2H^^{2vc\d + Svcvd) 



[Hab - 9abH f 



'9AB 



P + l 



2 (H^c + 2k)K - TJ^^cAp - ^H^''{r^K\c)\D - ^KK 



^[r^ikK)icf^ + ^Kick^^ + 2H^''^[iHFCv^ - cvf)9de + vfHde] + ^Hfc^H^^'d 



P + 1 <^'i? 



9DF I ^2 



— ) + 2{2vd\f + Svdvf) 



I I AB 



2i 1"- ' 

^-^1 \ HABh^ \c + ch(A\B) — 2H'~" i^AhB)\c + 2{Hab\c - Hc{A\B))h' 



c 



(85) 
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9AB 



H 



CD 



{r'hc)\D - H^ch'^iD + '2{H'-'dic - H'^ c\d - K^dW 



We have employed that the sums in I and I are symmetric to simphfy the form of the sources. Although each individual 
source (<')5'p!"™ is not symmetric under the interchange (/, ffi) [l, in), their sum is symmetrized. We have also tried 
to simplify the expressions as much as possible by using the GS scalar 11 = t^^ {r~'^hA)\B- In addition, we have 
defined Habc = Hab\c + Hac\b - Hbc\a- 

On the other hand, the source of A[Ga6] can be decomposed as 



I I A 



+ ^ofA/'i 2^bc\aH"^ + H"^ {Hbc\a - Hab\c - Hbcva) + -^{H" b\c ~ "^H" c\b)Ha^ 



+ hBVA)h'' ~ h'^hMB)} 



{kk)\A + T^K^aH^b 



Hh[A\B] + hBVA)h^ - h^hiA\B) + r^hAir-^h^)\ 



I I A 



-E- 



E^%i[HABh'' + h^HAB] 
P + l 



Olfh 
1 irhl 



-HABh" + h^'HAB + [H'^B - 2K)hA - hAin'^B + 2K) 



2 , 



A-t^lB 



(86) 



(87) 



2r-^{r^h^\Bk\A - 2r-^h^{r^K\B)iA + -hA 2^ [<^>/if ] + K{r^ '^'i? - - 1) + (r^A')!^^ 



2r\'''^IlHAB\c - r^ABUiHc^^'' - 2i?^^|c.) + 2r-hAB{r^tl)^cH''^ 



Finally, the manipulations for A[Ga6] are more complicated, involving up to 3091 terms in some intermediate steps. 
The resulting expression can be organized in the following four sources: 



1 



E^ii^A H^'^Hab - t^H^'aH^b + 2h^\^hA\B - -{rh^)\A{rh''\ 



2h^h'' 



-E 



2RaB + iVAVB - 9AB 

i^ + iP + i 



P + l-l-2rr\^ 



e^^hA{r'li\B 



+ 2^ i2HAC\B ^ 3H 



h^hA + 2{H^A - km^^'lBC - H^B^^'c) 
, 1 



AB\c) -+2^211 B\A — H a\b){2H \c — H c ) 



H 



AB 



Hal 



P + l 



- 2vbH''c\a + mAcv'')\B + -{rHAcy'^VB - 2Hab\cv^ 



{ — )Qmmm 

I I 



2[H^Br-\r^K)\B]\A + A\Br-\r'' K)^'' + '^^RH^aK - AH^'^Kvavb - K^^K^, 
^E^^i{ (g^^'H^c H^"") hA\B + (k\B + ^H^[A\B] H^aVb) hA, 



1 



tAB 



p + l -I 



h^hA\+El\t^H^^HA 



f ?n / 



I I 



-2iEl^i{kh^)\A - ^Ef^i [2{H^BhA)\B - H^MBh""] 



(88) 



(89) 

(90) 
(91) 



In spite of the obvious increase of complexity from the 
first to the second-order equations, we want to stress that 
the final expressions given above are still manageable and 
fully general, except for the choice of RW gauge. They 
can be particularized to the case of any spherical back- 



ground, dynamical or not, containing any type of mat- 
ter, and expressed in any kind of background coordinates 
(polar-radial, null, comoving, etc.) 

In situations with just a single first-order perturbation 
we will have / = I and rh = m. In these circumstances. 
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the iJ-coefiicients vanish for odd I. This implies that 
I + 1 — I is always even and therefore the sources '■^^S are 
never excited in the polar equations (|7S|l ~ (|7^ . neither are 
the sources '"'"'S' in the axial equations (|80|l and 181|) . In 
particular, '-"^Sab and '"'S* are never excited. This has 
been the case encountered in many previous investiga- 
tions in second-order perturbation theory. For instance, 
for the perturbations of a slowly rotating star a single 
axial / = 1 mode has been assumed and the studies 
of the close-limit black hole collisions have considered a 
single polar I = 2 first-order perturbative mode 0|. It 
will be very interesting to discuss what kind of interac- 
tions between modes are excited through the other types 
of sources, an issue that we plan to analyze in future 
investigations. 

D. Energy- momentum conservation 

A complete set of evolution equations is obtained only 
after specifying the particular type of matter content of 
the system (including as such the vacuum) . Some simple 
systems like scalar fields of perfect fluids are completely 
defined dynamically by energy-momentum conservation, 
but this is not the case in general. However, we can gen- 
erally analyze the consequences of perturbing the matter 
conservation equations, as well as use this analysis as a 
check of the perturbed Bianchi identities, and hence as 
a consistency check of the sources given in the previous 
subsection. 

Like every object in perturbation theory, the energy- 
momentum conservation law can be expanded into a hi- 
erarchy of linear equations, all sharing the principal part, 

A"[r^,-^] ==0. (92) 



At zeroth order we use the decomposition (|18|l provided 
by the background, which leads to a nontrivial relation 

QvA = ^{r^tAB)^'' ■ (93) 

Given its vectorial character, the energy-momentum con- 
servation equation can be decomposed into three geomet- 
ric parts at higher orders: a vector equation in the polar 
sector and two scalar (one polar and one axial) equations. 
At first order, those equations can be written in compact 
notation as 

LAnr,'''hr]=0, (94) 

L[f^vr/^*/in = o, (95) 

L[^^Vr. "V] = 0. (96) 

The operators La, L, and L are defined in Appendix IdI 
The second-order perturbation adopts the same form 
with additional quadratic sources: 

L^['=vr/'V] + ^ ^ '^'xpf rA = 0, (97) 

Li ^n,m 

L[f=vr, '^-'hr] + E E '^'^ff r = o, (98) 

l i m,m 

L[f=vr, '^'hr] + E E '-"^rfT = o. (99) 

l i rn.m 

Such sources can be computed starting from Eq. H92|l . 
decomposing it using formulas H6U|I and (|61|) . and finally 
applying the tools that we have developed to deal with 
products of harmonics. The result, with the same nota- 
tion used for the sources of the main equations, is 



{+)-j-rnmni 



I i I A 



-^E^%}^2r^''ir-^hB)\A - A - '^Q{r~^hBh'')\A + 2{^^Ah'')\B 

rHA''{r-'h^hc)\B - 2Ch^h^tAB)\c ~ rhBc{r~'h^h^')iA 



(100) 



•Ohn 
[fhl 



l^^^y^A-^A + 2*A^A'|B + 2QrK{r~^k)^A - 2KK^^tAB - 2r^^{r~^K)^A 



2H''^tABK\c - ^^^Hbc\a + ^A^H^c\B - ^{r^i'ABH''^) 



\c 



H 



BC 



4:Hc^\{AtD)B - H^D\ctAB - HBC\DtA^ + —{r'^HBDtA^)\C 



{ — yrmmm 



{+y-r7n7nm 



,2^1 iml 



H^Bi^A + {2H^^hc - H^ch^)tAB - 2/i^*AS 



(101) 



2r^^''{r-^hB)iA-'2{^Ah'')\B \, 



^E-^\'^il^2hAi''' + 2i,Ah^ + 2(^/i^)|A - ifH^A - Qh^hA - 2h^h''tAB 
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{ — y-rmmm 



2(r^^'Ai?^^)|B + 2r 



tpAh'^ " ^Ah"^ - h'^h^tAB + ^h'^hA 



+ ([- 1)([+2)A'^' 



■^'^K\A --^K- + QKK + ^^^H'^A] 



2HBct^A - ^AB - ^HaB + 9AB - 



-^2 



(102) 



2„7;Ar>. I j^/ri^27;AN. . , ^r^„2iTABr^_ , „,2 £tA , _ f tb Q tb 



2r>^X|^ + if(Qr^/i-^)|^ + {Qr'H^^hB)\A + r^i^^A|B ( - f 



r 



VI. CONCLUSIONS 

Perturbation theory has been highly successful in Gen- 
eral Relativity and nowadays still plays a relevant role in 
the simulations of the dynamics of many physical and 
astrophysical systems. On the one hand, it allows us to 
interpret certain dynamical processes as the evolution of 
a perturbative mode of a simpler process, or as the in- 
teraction among several of such modes. On the other 
hand, using perturbation theory we can evolve systems 
in which very different physical phenomena are happen- 
ing simultaneously with so distinct space/time scales or 
amplitudes that a numerical simulation would fail to fol- 
low them all with a satisfactory precision. Even so, one 
must be aware of the problem of linearization stability 
when perturbing a given spacetime in General Relativity, 
namely, it may happen that a perturbatively constructed 
spacetime with parameter e does not correspond to an ex- 
act family of solutions g^u{^) of the Einstein equations. 
This problem becomes more important when perturbing 
highly symmetric backgrounds |4ll |. 

However, going beyond first-order perturbation theory 
has not been feasible until very recently except in very 
specific situations. The present work proposes a system- 
atic approach to high-order perturbation theory in Gen- 
eral Relativity, based on the combination of a good choice 
of the theoretical formalism employed for the descrip- 
tion of the problem (implementing symmetry reductions, 
covariant notation, and other nice features) and the in- 
tensive use of abstract computer algebra to manipulate 
the enormous expressions that unavoidably appear in this 
field. 

We have first given a number of formulas which permit 
to compute very efficiently and at any order the pertur- 
bation of all relevant curvature tensors in General Rela- 
tivity, and have implemented them in the Mathematica 
package xPert. These formulas could be used in very dif- 
ferent areas of gravitational physics, including theories 
that depart from standard General Relativity (like in the 
case of models with extra dimensions, curvature correc- 
tions, or in braneworld scenarios). 



We have then generalized to higher orders the well 
known GS formalism for nonspherical first-order pertur- 
bations of a spherical spacetime. This formalism is con- 
sidered to be optimal for the perturbative study of a num- 
ber of astrophysical scenarios of interest, except for the 
possible construction of master scalars and equations re- 
ferred to in the Introduction. The generalization of the 
GS formalism put forward here will make the perturba- 
tive analysis even more powerful, leading to more precise 
results and allowing to describe interactions between dif- 
ferent modes. 

With this purpose, we have constructed a generaliza- 
tion of the GS harmonics, which turn out to be very 
closely related to the Wigner rotation matrices (spin- 
weighted harmonics in the General Relativity commu- 
nity). We have also obtained a general formula for the 
product of any pair of them, and implemented all these 
results in another package called Harmonics. 

In addition, we have computed all the equations of 
the generalized GS formalism at second order (includ- 
ing those of energy-momentum conservation) and sim- 
plified them to a form manageable enough as to allow to 
write them down in this paper. These equations are com- 
pletely general except for the restriction to a spherical 
background: they can be used with any background, dy- 
namical or not, they can be coupled to any matter model, 
and they have been given in covariant form, so that any 
coordinate system can be used on the background mani- 
fold. We have presented the results in RW gauge, which 
is well posed at all perturbative orders. The original GS 
formalism also employs gauge- invariant variables, and we 
arc currently working in this direction, with the aim at 
determining the gauge-invariants at second order in full 
generality. 

The second-order equations are essentially the same 
as the first-order equations, but they also include com- 
plicated quadratic sources. We have disentangled the 
structure of these sources and shown that, in previous 
investigations considering just a single first-order pertur- 
bative mode, many of such sources were not excited. We 
plan to study the role of these new sources in future work. 
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Let us conclude remarking that, once the gauge- invari- 
ants be determined, the formahsm and our computer- 
algebra tools implementing it will be ready to be ap- 
plied to problems of astrophysical or conceptual rele- 
vance, such as the emission of gravitational radiation in 
the collapse of a rotating star, or the coupling of pertur- 
bativc modes of a critical spacctimc. 
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APPENDIX A: SPHERICAL FUNCTIONS 

Several conventions arc employed in the literature for 
the special functions used in the theory of representations 
of the 3-dimensional rotation group. Here we mostly fol- 
low the conventions of Edmonds (3^ . 

The spherical harmonics Y™{9,(p) are 



Y, 



I I")' 



I {21 + ~ my. 
Att (l + my. 



P,"(cos6')e""'^, (Al) 



where P™ is the associated Legendre function 



2H\ 



The Mathematica functions SphericalHarmonicY and 
LegendreP are indeed those defined above. 

For a rotation of the reference frame described by the 
Euler angles {a, (3, 7), the convention adopted for the uni- 
tary matrix in a representation of SU{2) IS 

'D'-IM^P.i) = e'"'"dS™(/?)e™''. (A3) 
The /3-transformation is given by 



(A4) 



E 



(-1)'-™ -"^Jil + my.il - my.{l + m'y.il - m')! 
[I ^ m' — (t)! (I ~ m — (t)! (to' -|- m -|- cr)! cr! 



xlsm-j (cos- 



where the sum ranges over those integers a for which the 
arguments of the factorials are all nonnegative. 



APPENDIX B: SYMMETRIC TRACE-FREE 
TENSORS 

Given any tensor Pii...i, over a vector space of dimen- 
sion d with a metric gij , we construct its STF part as 



iSTF 



[Pii...iiJ — 

[;/2] 

E(m) ^ j-^ 



(Bl) 



with Si^,,,i^ = T(^ii...i,) and [1/2] the integer part of 1/2. 
The coefficients of the expansion are determined by the 
trace-free condition, and are given by 



'■Ld 



II 



r[l + d/2-l- m] 



(-4)™to!(; - 2m)! T[l + d/2-l] 



(B2) 



In our case, d ~ 2 for the unit sphere. These formulas 
allow us to compute any of the ZPai...as in terms of the 
derivatives yj™:ai...as or viceversa. Note that derivatives 
of Y™ with indices sorted differently are not equal, but 
can be transformed into a term with the desired order of 
indices plus terms with a lower number of derivatives. 



APPENDIX C: PURE-ORBITAL HARMONICS 

In this Appendix, we briefly discuss the formalism of 
pure-orbital tensor harmonics and its relation with the 
pure-spin harmonics. Firstly, one constructs pure-orbital 
vector harmonics by composing scalar harmonics of an- 
gular momentum I with a set of 3-dimensional vectors t"^i 
which transform under a representation of spin 1 |37j| : 



or., 



+1 
E 



f^m — 7n m m \rni — m 4_m 
1 J 



(Gl) 



with j = ^ — 1, Z, Z -f- 1 and |to| < j. The vectors f^i are 
defined in terms of a fixed orthonormal Gartesian basis: 



t 



±1 



V2 



(G2) 



(The index i is an abstract index on the manifold 
with Euclidean metric, in which is embedded.) These 
vector harmonics Of"^ i transform under a representation 
of total angular momentum j and their Gartesian com- 
ponents arc eigenfunctions, with eigenvalue 1(1 + 1), of 
the S'^-Laplacian (also called orbital angular momentum 

M) 



L 



2 — „2f72 



r^V^ + 9,(r^a,) = -7'^''VaV6 



(G3) 



Pure-orbital vector harmonics are however not trans- 
verse to the radial direction. 



y^OOO \rn-l 

'^3 II ^3 ■ 



(G4) 
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Therefore, one must take certain linear combinations can- 
celing their radial contribution to get the GS harmonics. 

Pure-orbital bases for higher-rank tensors can be con- 
structed recursively from the above vector basis. The ba- 
sis for STF tensors with s indices and well-defined spin 
s can be built by composition of the bases with s' and 
s — s' indices (with any < s' < s) as follows: 



which is a sum with only two contributions, /' = / ± 1. 

Pure-spin harmonic tensors, with indices on S^, can be 
obtained as linear combinations of the pure-orbital ones: 



j+s 



y;"'a,...a.^ J2 V2m4,,Of™a,...a.. (CIO) 
l=j-s 



s — s' s i^,. 



(C5) 



From this, we construct orbital harmonics with s indices: 

+s 



'-^l 41... 4s 



'm—m m m v^m— m j.m 
/ s j ^ ii...isi 



(C6) 



which are normalized so that 

dn (pf"\,...^Sof"'\,...,^ ^ <5ii'%'<5w. (C7) 

The symbol * denotes complex conjugation. 

One can obtain the following multiplication rule by 
using Eq. (|38|) and the formulas available in the literature 
for the composition of three angular momenta [s^ : 



(2/-f l)(2Z' + l)(2j + l) 

I'l" 



An 



I'+l 

E 

l" = \V-l 

j'=i"-i 



(C8) 



rn-\-m 



where W is the Racah coefficient [Sg . From this formula 
it is possible to compute the product of any two orbital 
harmonics using the Leibnitz rule, but here we prefer 
to show their connection with the pure-spin harmonics, 
which obey a much simpler multiplication formula. 
The radial component of the orbital harmonics is 



^(2s + l)(2/ + l) 

Ecf??,M/(.-i,i,j,/;.,norn...4.-i. 



(C9) 



Asking them to be orthogonal to the radial direction, we 
find from Eq. (|C9(I that the coefficients of this expansion 
obey the "second-order" recursion relation 



4!^ 



I is-i+m+j + i-s) 



(Cll) 
(C12) 



^ is + l-j + l){l + j + 2 + s)- 
Since the odd-^ coefficients decouple from the even-Z ones, 
there are two independent solutions. Fixing a normaliza- 
tion constant and employing definition (|47|) . we get 







4tt 



V2ITT V 2j + 1 



(C13) 



±s,m 



which define, respectively, the harmonics y, 

Finally, using formulas l)C6|l and H42|l , and the relation 
(I43|l (with s ~ 0) between rotation matrices and spherical 
harmonics, one can express the pure-spin harmonics as 



y. 



(CM) 



E 



(0,0,0) 



APPENDIX D: GS OPERATORS 

For the sake of completeness, we include in this Ap- 
pendix the metric part of the GS equations at first per- 
turbative order and in RW gauge. This is also the linear 
part of the perturbation equations at any order. As it is 
usually done at first order, we obviate here the labels /, 
m, and n of the perturbations on M^. 



(l-l)(l + 2) r r 

{vbK\a + vaK\b + K\ab) + 9AB 



Had + vc{H^B\A + a\b ~ Hai 



r-3(^3^|C)|^_ 



G-i)a + 2) 

2r2 



K 



2r2 

{H^'bva - H^'bia + Ha'^ib - K\a), 



(Dl) 



(D2) 



15 



ITT T^„ \(i)r>^B l{l + 1) . ttA B 

[Hab ~ KgAB) K — H a + n a\b 



2H^B\Av'' + H^A\Bv'' - H^''\AB + K\^a + 2K\AV^'^. 



G -!)(/ + 2) 1 



2r2 



hA 



he 



1^. 



(D3) 

(D4) 

(D5) 
(D6) 



On the other hand, the hnear parts of the equations for the perturbations of the energy-momentum conservation 
law, also in RW gauge, are the following: 



^l±ll^^-2vA^ + ^{rHABr 



1 ,BC TT ''^ 

2^ J^BC\A - 



■Qir-^K)\A - l:tABH^'c^'' ~ tABK^"" + \(rHABll''^) 



1 



a -1)^ + 2), 



Q 1 



I 



(D7) 

(D8) 
(D9) 



APPENDIX E: COMPUTER IMPLEMENTATION 

Three Mathematica packages have been constructed 
during the course of this investigation, respectively called 
xPert, Harmonics and xPertGS. This Appendix describes 
what is included in them. The three packages are based 
on the package xTensor |42l | for abstract tensor compu- 
tations, written by one of us, and also freely available. 

The package xPert implements the equations of Sec. ^1 
In particular, there is a command Perturbation which 
plays the role of A and is able to give nth-order per- 
turbations from those at order n — 1. The command 
Perturbation [ex^r, ri] is recursively computed via 
Perturbation [Perturbation [eccpr, n-l'i , i] . There 
is then the command GeneralPerturbation [expr, n] , 
which implements the expansion formulas for the metric 
(Q, the inverse metric (jS)), the Christoffel symbols (|12|l . 
the Riemann tensor H15|l . and the Ricci tensor and scalar. 
The general expansions are much faster than the recur- 
sive procedures, and for instance it is possible to produce 
within seconds the expressions of fifth-order perturbation 
theory with a small PC. 

The package Harmonics implements all the structures 
that have been defined on in Sec. IIIII Sec. IIVI and 
Appendices IXI IbI and lUl of this paper. The commands 
PureOrbital [j', I, m][a, h, . . .] and PureSpin , 
±1, m\ [a, b, . . .] give, respectively, any pure-orbital 
or pure-spin harmonics, both in abstract form (in terms 
of the bases t or m) or providing their components in any 



coordinate or noncoordinate basis. The generalized GS 
harmonics Z and X have also been defined, incorporating 
all their symmetries and properties (|31|l - l|37|) . as well as 
the product formula lfCT|) . 

The package xPcrtGS needs the previous two packages 
and has three sections, closely following the mathemati- 
cal structure of this article: 

1. Spherical background. The reduced manifold Af^ is 
defined with its metric g [-A , -B] (and correspond- 
ing derivative CD [-A] ). The scalar field r [] is also 
defined. The manifold is constructed as a prod- 
uct of and 5^; all tensors on can be block- 
decomposed in their respective parts on those sub- 
manifolds. 

2. GS first order. First-order perturbations are de- 
fined and their equations computed using the gauge 
invariants (the total execution time is of the order 
of two minutes in a PC). 

3. GS second order. Second-order perturbations are 
defined and their equations computed and manip- 
ulated in order to simplify them as much as possi- 
ble (the total execution time is of the order of two 
hours). 

The packages xPert and Harmonics can be freely down- 
loaded, under the GNU General Public License, from 
|http: //metric . iem. csic . es/Martiii-Garcia/xAct/xPert| 
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